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by 


Charles  E.  Blair 
Department  of  Business  Administration 
University  of  Illinois 


Abstract 


We  consider  methods  of  locating  a  completely  labelled  simplex  of  a 
triangulated  unit  simplex.  We  show  that  a  class  of  algorithms  which  includes 
the  usual  pivoting  methods  may  be  forced  to  evaluate  a  positive  fraction  of  all 
the  vertices  of  the  triangulation.  An  algorithm  which  does  not  have  this 
property  is  briefly  described. 


Let  T  be  a  triangulation  of  the  unit  d-dimensional  simplex.   If  L 
is  a  proper  labelling  of  the  vertices  of  T,  Sperner's  lemma  guarantees 
the  existence  of  a  completely  labelled  d-dimenslonal  simplex.   It  is 
well  known  that  the  problem  of  locating  a  completely  labelled  simplex 
is  essentially  equivalent  to  the  computation  of  an  approximate  fixed 
point. 

We  consider  algorithms  A  for  locating  a  completely  labelled  simplex 

of  T.   By  an  algorithm  we  mean  a  deterministic  procedure  which  asks  for 

the  labels  of  various  of  the  vertices  of  T,  in  such  a  way  that  the  next 

vertex  for  which  a  label  is  demanded  is  determined  by  the  labels  already 

known.   Eventually  A  will  locate  a  completely  labelled  simplex  and  halt. 

Let  N(A,T,L)  be  the  number  of  vertices  whose  labels  are  demanded  by  A 

when  used  on  (T,L).  We  are  interested  in  W(A,T)  =  max  N(A,T,L)— the 

L 
most  labels  A  can  be  forced  to  evaluate  by  a  malicious  labelling  L. 

In  the  terminology  of  [3],  we  are  regarding  the  labels  of  T  as  given 
by  an  oracle  and  W(A,T)  is  related  to  the  oracle  complexity  of  A.  Worst- 
case  analysis  of  other  problems  have  been  treated  in  [4],  [5],  and  [6]. 

Most  of  the  popular  methods  of  finding  completely  labelled  simpllces 
proceed  by  locating  a  completely  labelled  (d-1) -dimensional  simplex  of  T 
on  the  boundary  B  of  the  simplex  and  then  moving  through  a  sequence  of 
almost  completely  labelled  d-dimensional  simpllces  of  T.   Our  definition 
of  a  local  algorithm  is  a  slight  generalization. 

Formally,  a  local  algorithm  is  one  that  generates  a  sequence  of  (d-1) 

dimensional  simpllces  S-,...S  of  T  and  vertices  ?^,...'P     such  that 

In  in 
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(i)    (S,,  P  )  is  a  d-dim  simplex  of  T. 

(ii)  S^^^   is  a  face  of  (S^,  p^)  and  p^_^^  i    (S^,  p^) 

(iii)  If  J  is  a  (d-1) -dimensional  simplex  in  B  there  is  a  labelling 
of  B  which  will  make  S^  =  J. 

(iv)  If  J  is  a  (d-l)-dim  simplex  of  (S^,  p^)  J  ^  S^  there  is 

a  labelling  of  p  which  will  make  S^_^^  =  J  and  (S^,  p^)  will 
not  be  completely  labelled. 

The  crucial  property  is  (iv) ,  which  essentially  says  that  the  next 
simplex  (S .  , ,  ,  P^.-i)  visited  by  A  is  determined  by  the  labelling  of  p^.  The 
algorithms  we  have  in  mind  are  those  in  which  S .  is  a  completely  labelled 
(d-1) -dimensional  face  of  (S^,  p^)  and  S^_^^   is  the  uniquely  determined 
other  completely  labelled  face.   It  is  clear  that  such  algorithms  satisfy  (iv) . 

Let  J„(N)  be  the  Whitney- Freudenthal  triangulation  of  the  d-dimensional 
unit  simplex  (we  consider  d  fixed  throughout).   See  [7,  p.  32-35]  for  a 
precise  description. 

Theorem  1;   There  is  a  K  such  that,  if  A  is  any  local  algorithm 

W(A,J2(N))  >.KN^   . 

Proof:   Define  a  snake*  of  length  M  as  a  sequence  Q^^,  Q2,...Qj^  of 
d-dimensional  simplices  such  that  (a)  Q.  shares  a  (d-l)-dimensional  face 
with  Q.., .   (b)  Q,  has  exactly  one  face  on  B.   (c)  the  unique  point  of 
Q.  .  not  in  Q.  is  not  in  Qi.'"Qi_i'  ^^^    the  only  vertices  of  Q^  in  B 
are  in  Q. . 

*This  term  is  taken  from  [2],  where  it  was  used  for  a  related  idea  in  a 
different  context. 
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We  claim  W(A,T)  _>   the  length  of   the  longest   snake   in  T   (this   is 
true  for  any   triangulation) .      By   (iii),    the  boundary  B  can  be  labelled 
so  as   to  make  A  choose  S^    as   Q^OB.      Then,    for   each   i,    p,    is   labelled 
so  as   to  make   S.^    =  Q.OQ.^.      Condition   (d)    insures   that   the  propriety 
of   the  labelling  does  not   interfere  and   condition   (c)    insures  p .    ?*  p , 
for   i  T^  j. 

To  complete   the  proof  we  must  construct  a  long   snake  in  J«(N) . 

Following   17],   we   identify  J^(N)   with   the  corresponding   triangulation  of 

"UjCN)   =  {(x,,...x,)    In  >  X,    >    ...   X,   >   0}.     We  construct  a  path  P   in  U,(N) 
a  id'      —     1  —  a  —  a 

and  obtain  the  snake  by   taking  Q^ ,   Q-,...    to  be   those   simplices  of  ^2^^^ 
such   that  P  passes   through   their   interiors.      Intuitively,    P  will  move 
through  the  layers  of  U,(N)    (N  1  X.    >^  N-2,    N-2  ^  X^  _>  N-A,...)    so   that 
P  passes   through  a  positive  fraction  of  all   the  simplices  of  J„(N). 

The  construction  of  P   takes   several   steps.      Let  0  <   f^    <   f-  <    ••• 
f ,   <   f ,   .    <  1.     For   each  1  ^  j  £  d  and  M  B    (2j   +  1)   mod  4  we  construct 
a  sequence  of  points  V(j,M)    in  U.(M  +  1)   whose  first  member   is 
(M  +  f .,  M-2  +  f .   ,,...M-2i  +  2  +  f, )    and  whose   last  member   is 
(2j-l  +  f  ..,,    2j-3  +  f.,...l   +  f „) .      Note   that  V(j,M)    is  an  ordered   subset 

j+1) .     We  will  denote  by  -V(j,M)    the  same  set   in  reverse  order. 

For   j=l  V(j,M)    is  defined  as  <M  +  f j^,   1  +  f^>'.   For   j   >   1  V(j,M)    is 
defined  as   (M  +  f .,   V(j-1,  M-2));    (M  +  f.   -2,   -V(j-1,  M-6)); 
(M+  fj  -  4,   V(j-1,  K-6));...(2j  +  1  +  f^,  V(j-1,    j-1)); 
C2J-1  +  f^^^,    2j-3  +  f^,...l  +  f2). 

Here   (x,   V(j-1,   y) )   denotes   the   sequence  of   points   in  U.(>H-1)  whose  Kth 
term  has  first  co-ordinate  x  and   2nd   through   j-th  coordinates  determined  by 
the  Kth  term  of  V(j-1,   y)),    and   the  sequences  are   concatenated. 
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Example;     For  f^  =    .li,   V(3,ll)    is   (11.3,    9.2,    7.1); 
(11.3,   9.2,   1.2);    (11.3,    7.2,   1.2);    (11.3,    7.2,    3.1);    (11.3,    5.2,    3.1); 
(11.3,   5.2,   1.2);    (11.3,    3.3,   1.2);    (9.3,    3.3,    1.2);    (9.3,    5.2,   1.2); 
(9.3,    5.2,    3.1);    (7.3,    5.2,    3.1);    (7.3,    5.2,   1.2);    (7.3,   3.3,   1.2); 
C5.4,   3.3,   1.2). 

We  can  now  construct  P.      Let  N-5  £  M  _<  N-2,   M  =   2d  +  1    (mod  4)  . 

P  is   the  polygonal   path  obtained  by  connecting   the  point   (N,   M-2+f        , 

M-4+f,    „,...M-2,+2+f,)    to    the  point    (M+f ,, . . .M-2d+2+f   )   and   then  connecting 
a— 2  d         1  a  1 

the  points  of  V(d,M)    in  follow-the-dots   fashion. 

P  generates  a  snake  of  length  >_  (1/2)  .^_j^v  ^   N     so   theorem  1   is 

proved  with  K  >_  (1/2)^'^^    ,^]^^.  ^      .      Q.E.D. 

A  similar  result  can  be  proved   for   the  Kuhn  triangulation  K2(N) . 
I  conjecture  that  for  any  local  algorithm  A  and  any   triangulation  T 
of  mesh  <_  f  W(A,T)    =  0(f'*^)  . 

However,    theorem  1  does  not   extend    to  non-local  algorithms. 

Theorem  2:     There  exist  non-local  A  such  that,   for  N=2   ,  W(A,J2(N))  1  KN 
Proof:      A  is  based  on  an  extension  of   Sperner's  lemma    [1,    theorem  4]: 
Let   the  triangulated   simplex  T  have   each  of   its  vertices   labelled 
with  {0,1,... d}   in  an  arbitrary  manner    (i.e.,   not  necessarily  a  proper 
labelling) .      Then   the  parity  of   the  nvmber  of   completely  labelled  d-dlmen- 
sional   simplices   is  equal   to   the  parity  of   the  number  of   (d-1)    dimensional 
slmplices   in  the  boundary  B  which  are  labelled  with   {l,...d}. 

The  propriety  of   the  labelling   insures   that   the  unit  simplex  contains 
an  odd  number  of    (d-1) -dimensional   simplices  labelled   {l,...d}   in  its 
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boundary.  A  divides  the  unit  simplex  into  two  simplices  W^ ,  W„ 
whose  interiors  are  disjoint,  and  such  that  the  hyperplane  W.^_  does 
not  pass  through  the  interior  of  any  d-dimensional  simplex  of  J^(N)  . 
Either  W.  or  W„  must  contain  an  odd  number  of  (d-1) -simplices  labelled 
{l,...d}  in  its  boundary.   The  extended  Sperner  lemma  then  implies 
that  a  completely  labelled  simplex  is  in  W^  in  one  case  and  W»  in  the 
other.  Then  A  divides  W.  (i=l  or  2)  into  two  simplices  W„  or  W,  in 
the  same  manner.  Again  the  extended  Sperner  lemma  tells  which  of  W_  or  W, 
must  contain  a  completely  labelled  simplex,  and  so  on.  At  each  step,  we 
have  a  simplex  which  is  known  to  contain  a  completely  labelled  simplex. 
This  simplex  is  divided  into  two  sub-simplices  and  the  extended  Sperner 
lemma  enables  one  to  determine  which  one  must  contain  a  completely  labelled 
simplex.  This  determination  requires  evaluation  of  the  labels  on  the 
boundaries  of  the  two  subsimplices  but  not  in  their  interiors.   After  dL 
bifurcations  a  completely  labelled  simplex  of  J„(N)  is  located. 

The  existence  of  hyperplanes  at  each  stage  which  divide  the  simplex 
in  two  and  do  not  pass  through  the  interior  of  any  simplex  of  the  tri- 
angulation  is  a  special  property  of  J»(2  ).   As  before,  we  regard  J„(N) 
as  a  triangulation  of  U,(N).  At  the  first  step  we  divide  U,(N)  using 
the  hyperplane  X.  +  X ,  =  N.   This  hyperplane  does  not  pass  through  the 
interior  of  any  simplex  of  J„(N).  Further,  the  affine  mapping  (X^, .  .  .x,)->- 
CN-X,,  N-X,  -  ,...N-X2,  N-X  )  carries  the  half  X^  +  X^  _<  N  onto  the 
half  X-  +  X ,  2L  N  in  such  a  way  that  simplices  of  J^CN)  are  mapped  to 
simplicies  of  J„(N).   This  means  the  two  halves  are  equivalent.   After 
K  <  d  bifurcations  one  obtains  a  simplex  equivalent  to  Uj(N)^X^  +  Xj_,  ..  <_  N. 
The  hyperplane  X^  +  X ,  ,  =  N  divides  this  simplex  into  two  halves  that  are 


-6- 

equivalent  with  respect  to  the  mapping  (}L  ,.  .  .X,)->-(N-X^_^,  N-X^_2 , . . . N-X^^ , 

3Lj.,...X,).   After  d  bifurcations  a  simplex  equivalent  to  U,(1/2N)  is 

obtained.   This  is  divided  in  half  by  X^  +  X ,  =  y,  and  so  forth. 

Each  face  of  J„(N)  contains  (  ,  ,   )  vertices.   Each  of  the  next  d 
I  d-l 

divisions  requires  computing  the  labels  of  the  vertices  on  the  dividing 

hyperplane,  which  is  <(  ,  ~  )  vertices.   The  problem  is  then  reduced  to 

N 
J.(y)  so  the  total  number  of  vertices  for  which  labels  are  demanded  is 

<  (2d+l)  C?*"?)  (l+(l/2)'^"-'-+  (1/2) ^'^"^  +  ...  This  establishes 
—        a— 1 

W(A,J2(N))  <_   (2d+l)  (l/l-(l/2)'''"^)  (^^"^  )  =  0(n'^"^).         Q.E.D. 

Concluding  Remarks : 

1.  A  more  detailed  analysis  shows  WCA.J-CN))  <_  3  (   ~  )  for  N  large. 

2.  For  N  not  a  power  of  2  J„(N)  can  be  embedded  in  J2(M)  for  M  the  next 
greater  power  of  2  and  a  bifurcation  algorithm  can  be  performed  on  J2(M). 
This  will  still  require  evaluation  of  0(M   )  labels. 

3.  The  function  W  is  an  artificial  measure  of  the  amount  of  time  required 

by  an  algorithm.   It  ignores  any  computational  effort  other  than  evaluations 
of  labels,  so  it  is  only  relevant  in  situations  where  individual  evaluations 
of  the  function  are  time-consuming.   Even  in  this  context,  the  complete 
bifurcation  algorithm,  of  Theorem  2  is  probably  impractical  for  d  >  3, 
but  some  hybrid  of  bifurcation  and  the  more  familiar  local  algorithms 
may  be  possible.   Further  discussion  may  be  found  in- [1],  which  discusses 
bifurcation  in  a  non-simplical  setting. 
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